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ABSTRACT 

Well-determined physical parameters of 130 W UMa systems have been collected from 
the literature. Based on these data, the evolutionary status and dynamical evolution 
of W UMa systems are investigated. It is found that there is no evolutionary difference 
between W- and A-type systems in Af — J diagram which is consistent with the results 
derived from the analysis of observed spectral type, M — R and M — L diagrams of W 
UMa systems. M — R and M — L diagrams of W- and A-type systems indicate that a 
large amount of energy should be transferred from the more massive to the less massive 
component so that they are not in thermal equilibrium and undergo thermal relaxation 
oscillation (TRO). Meanwhile, the distribution of angular momentum, together with 
the distribution of mass ratio, suggests that the mass ratio of the observed W UMa 
systems is decreased with the decrease of their total mass. This could be the result 
of the dynamical evolution of W UMa systems which suffer angular momentum loss 
(AML) and mass loss due to magnetic stellar wind (MSW). Consequently, the tidal 
instability forces these systems towards the lower q values and finally to fast rotating 
single stars. 

Key vifords: binaries: close - stars: magnetic fields - stars: mass-loss- stars: evolution 



1 INTRODUCTION 

, W UMa stars are short-period, dumbbell-shaped binaries in 
■ which both stars are in contact or overflowing their Roche 
limiting surface. They are classified in to A- and W-type sys- 
tems on the base of the light curve (|Binnendiiklll970D . A- 
type systems are showing primary minima due to the eclipse 
of the larger, more massive component while the reversal is 
true for W-type systems. However, a satisfactory theory for 
the origin, structure and evolution of the W UMa binaries 
has not yet been suggested. The traditional view for the ori- 
gin of contact binaries is that W UMa systems are formed 
from detached binarie s of comparable periods through or- 
bital decay by AML (IVilhulll98ll: iRahunen fc Vilhulll98ll : 
iPaczvhski et al" I I2OO6I ). IPaczvhski et all (120061 ) have found 
too many contact binaries in comparison with the possible 
detached progenitors having orbital periods about one day. 
Then, W UMa systems seem to appear out of "nowhere". 
Even if it is possible that W UMa systems originate from 
detached progenitors through AML due to MSW if they 
are formed from detached binaries with periods P < 2.24 
d (corresponding to a maximum lifetime of about r < 3.23 
Gyr in pr e-contact stag e) and the lifetime of them are 
very long (|Li et al.|[20OT ). Otherwise, one should not find 
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so many contact binaries through ASAS. The excess of W 
UMa stars implies o ther channel of formation. Recently, 
IPribulla fc Rucinskil l|2006h have found that up to 59 per 
cent of W UMa bina ries ha v e com panions. This opens up 
a possibility that the iKozail (|l962l ) cycle operates in some 
triples. But, the relative importance of the traditional origin 
of W UMa binaries through AML, and through the Koza i 
cycle, will be required further study l|Paczvhski et al.ll2006l ). 
With regard to the lifetime of W UMa binaries, the different 
lifetimes have been derived for W UMa systems by various 
authors based on different m aterials, suc h as 0.1 — 1.0 Gy r 
jGuinan fc Bradstreetl Il988l ). 1.61 Gyr (iBilir et al.1 l2005l) 
4.7 Gvr (|van't Veerlll996l)T7.2 G yr l|Li. Han fc Zhangll2005l ) 
and >5.68 Gyr (|Li et al.1120071 ). 

Wit h regard to the evolutionary statu s of W UMa 
systems, iHilditch. King fc McFarland l|l988t) have argued 
that the primary components of the shallow-contact W- 
type systems are unevolved main-sequence stars while those 
of A-type ones are near to the terminal-age main se- 
quence. Meanwhile, there are other physical differences be- 
tween these two subclasses. A-type systems usually have 
an earlier spectral type, a larger mass, a higher luminos- 
ity and a smaller mass ratio. The degree of overcontact 
is larger and a thick common envelope is present. W- 
type systems are showing a later spectral type, a smaller 
mass, a lower luminosity, and a larger mass ratio. They 
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have a shallow convective envelope. So A- and W-type 
systems are usually assumed to be in the different evolu- 
tion schemes or evolutionary status. However, the observed 
spectral type does not indicate a difference in evolution- 
ary s tate between A- and W-type subgroups (IVan Hammd 
1 19821 . and references therein). In addition, some interest- 
ing ideas have been proposed concerning the evolutionary 
link between A- and W-type systems. Some theories sug- 
gest the possibility that W-t ype systems evolve into A- 
type through mass exchange (iHilditch. King fc McFarland 
1 19881 : Im aceron i. Milano & Russo"l98^ while ot hers suggest 
the opposite (Rucinski 1985; Gazoas fc Niarchos 2006). Here 
three questions arise: do both subtypes have the same ge- 
netic origin? regarding A types and W types, is one type 
the progenitor of the other one and what will be the final 
outcome of the evolution of each grou p? The above ques- 
tions are still open for investigations l|Van Hammd 1 19821 : 
iGazeas fc Niarchos! [20061 ). 

The components in a W UMa system have nearly equal 
surface temperatures in spite of their often greatly differ- 
ent masses. The canonical explanation of this phenomenon 
is that a large amount of energy is transferred from the 
more massive to the less massive component. However, the 
problem whether or not W UMa systems are in thermal 
equilibrium also remained an open question. It is interest- 
ing that there is controversy about it at present. Some au- 
thors claim that a contact configu r ation can achieve ther- 
mal equilibrium, such as IStepieiil (|2006l ) argued that W 
UMa systems can achieve thermal equilibrium if the sec- 
ondary, which was the former primary, is more advanced in 
its evolution, in analogy with the 'Algol paradox', so that 
W UMa systems should not undergo TRO. Ho wever, it is 
not c l ear that th i s asse rtion is generally correct l|Paczvhskil 
I2OO7I ) . iPaczvnskil (|2007l ) proposed that this is likely applica- 
ble to systems with extreme mass ratios. In addition, other 
authors claim that a contact configuration in thermal equi- 
librium is not po ssible so that W UMa systems are forced 
to un dergo TRO (iLi. Han fc Zhandl2004l : lYakut fc EggletonI 
l2005l : |Paczviiski et ahlbOO^ 

In this work, we have collected the well determined pa- 
rameters of 130 W UMa systems (69 W- and 61 A-types) 
from the literature. Based on these data, the evolutionary 
status and the dynamical evolution of these binary stars are 
discussed in the following Sections. 



2 STATUS OF W UMA SYSTEMS 

We collect the well-determined physical parameters of 
78 W UMa con t act bi naries from the c ompile s of 
Yakut fc EggletonI (120051'): lAwadalla fc Hannal (|2005l ) and 
Maceroni fc van't Veed l| 19961 ). In addition, we collect the 
new physical parameters for 52 W UMa systems from other 
literature (listed in Table 1). Based on these data, the evo- 
lutionary status of W UMa systems has been analyzed. The 
relations oi M — R and M — L for W- and A-type systems, 
together with the relations of M — R and M — L of zero age 
mam sequence (ZAMS) stars (|Tout et al.lll996l ). are plotted 
in Figures 1 and 2, respectively. It is seen in Figure la that 
most of the primaries of W-type systems are above ZAMS 
and seem to have evolved away from ZAMS, only the pri- 
maries of 13 W-type UMa systems (GZ And, 44i Boo, AC 



Boo, EF Boo, VW Cep, TW Get, TX Gnc, SW Lac, ET 
Leo, ER Ori, BB Peg, HT Vir, and BD-f42° 2782) are un- 
der ZAMS. However, as seen from Figure lb, apart from the 
primaries of 13 W UMa systems mentioned above, the pri- 
maries of 27 other W-type systems are found to have moved 
under ZAMS. Meanwhile, it is seen in Figure 2a that most of 
the primaries of A-type systems are above ZAMS and seem 
to have evolved away from ZAMS, only the primaries of 5 
A-type systems (V417 Aql, BI GVn, V899 Her, VZ Lib, and 
EQ Tau) are under ZAMS. But as seen from Figure 2b, apart 
from the primaries of 5 W UMa systems mentioned above, 
the primaries of 23 other A-type systems are found to have 
moved under ZAMS. This indicates that the evolutionary 
status of the primaries of the W UMa systems (including A- 
and W-types) predicted by M — _R relation seem to be more 
advanced than that predicted by M — L relation. Meanwhile, 
both AI — R and M — L diagrams show that the secondaries 
of A- and W-type systems appear to have been evolved away 
from ZAMS. The appearance might be as a result of the en- 
ergy transfer from the primary to the secondar y in W UMa 
systems. According to iHazlehurst et al.l (|l977l ). the effects 
of energy transfer on the radii and effective temperatures of 
the primaries of W UMa systems can be given by 



Alogii(t) = / ip^{t-T)EdT, 
Jo 

Alogre(t) = / i,^{t~T)EdT, 

Jo 



(1) 

(2) 



where E is the rate of energy transfer, tp^ and ■j/'t 
response functions for the radius and effective temperature 
with respect to energy transfer. The luminosity L oc R?T^, 
the effect of energy transfer on the luminosity can be given 
by 



AlogL(t) = 2Alogi?(t) +4AlogTc(t) 



(2i/)R-h4?/>^)(t-r)Sdr. 



(3) 



The primaries of W UMa systems always lose the energy 
in the common envelope, the energy loss would lead the 
effective temperature and radius of the primaries to be de- 
ceased. Accordin g Eqs. (1) and (2), %!}„ and Vt should be 
positive (also see IHazlehurst et al.lll977l ). Therefore, the re- 
sponse function (i/)^ = 2i/)j^ -I- ^ijj,^ ) for the luminosity with 
respect to energy transfer is larger than that (i/i^ ) for the 
radius, which implies that the luminosity of the primary is 
affected by energy transfer more greatly than its radius in 
the W UMa systems. Although the stellar evolution would 
lead to continuous increase of luminosity and radius, the en- 
ergy transfer (loss) would lead their luminosities and radii to 
be decreased and attempt to draw them back to ZAMS line. 
Since the luminosity of the primary is affected by energy 
transfer more greatly than its radius, the positions of the 
primaries should be drawn back by energy loss more greatly 
in M — L diagram than m M — R diagram. It is a main rea- 
son why L values of the primaries appear to be lower than 
ZAMS line while R values of the primaries are higher than 
those of the systems on ZAMS. This suggests that the sta- 
tus of the primaries of W UMa systems indicated by M — _R 
diagram seems to be more advanced than that predicted by 
M — L diagram. 
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Figure 1. Relations of mass-radius and mass-luminosity for W- 
type W UMa systems. Open and full dots represent the less mas- 
sive and more massiv e componen t s, res pectively. The solid lines 
represent ZAMS from lTout et al.l l[l99i). 
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Figure 2. Relations of mass-radius and mass-luminosity for A- 
type W UMa systems. Solid lines and symbols are the same as 
Fig. 1. 



Figures 1 and 2 also show that both radius and luminos- 
ity of the secondary of each W UMa system are evidently de- 
viated from those of ZAMS stars, and that the evolutionary 
status of the secondaries appears to be more advanced than 
that of the primaries. Most probably, this must be the result 
of the energy transfer from the primary to the secondary in 
the W UMa systems. In general, the energy transferred from 
the primary to the secondary in the envelope of W UMa sys- 
tems is severe or even tens of times more than the nuclear 
evolutionary effect on the secondaries so that the positions 
of the secondaries are significantly different from those of 
ZAMS stars. 



3 THE POSSIBLE OUTCOME OF W UMA 
SYSTEMS 

The angular momentu m distribution of 78 W UMa sy stems 
had been discussed by iMaceroni fc van 't Veeil { 1996f). It is 
found that most of A-type systems seem to have no evo- 
lutionary link with W-type ones because A- and W-type 
systems can be divided into two separate parts in M — J 
diagram by a line which gives the angular momentum of a 
system with just in contact main sequence components and 
a primary mass of 1.35 Mq. The orbital angular momentum 
of a binary can be written as 



(5) 



, Ml Ma .2 

■JoYh — A LOo 



(4) 



where M and Mi, 2 are the total mass and the masses of 
both components, respectively, A the orbital separation and 
Wo the orbital angular velocity. The spin angular momentum 
of the binary is 



where i?i,2 are the radii of the primary and secondary, fci,2 
the ratios of the gyration radius to the stellar radius for 
both stars and lus the spin angular velocity. In general, the 
rotational periods of the components are synchronized with 
the orbital periods for W UMa systems, due to strong tidal 
interaction (i.e. ujs = ^j^o = ij^)- If the ratios of the gyration 
radius to the stellar radius are assumed to be equal (i.e. 
fci = ^2 = fc^), the total angular momentum of the binary 



J = 



■ M1M2 



+k'{MiRl + M2RI)]u 



(6) 



>- M 

T aking = 0.06 as in iRasiol l| 19951 ) and iLi fc Zhan j 
(|2006t ). the angular momentum distribution of 130 W UMa 
systems with well-determined parameters is investigated. A 
very interesting feature is shown in Figure 3, where the to- 
tal angular momentum of 130 observed W UMa systems 
based on equation (6) is plotted vs the total mass with full 
and open dots for A-types and W-types, respectively. Fig- 
ure 3 confirms the fact that the total angular momentum 
(J) decreases by decreasing total total mass (M) of W UMa 
binaries. 

In order to obtain the information on the dynamical 
evolution of W UMa systems from the angular momentum 
distribution, it is necessary to give a theoretical distribution 
of angular momentum vs the total mass for W UMa systems. 
If u is in units of day~^, M in Mq and A in Rq, the Kepler's 
third law can be written as 



. = 54.23(|^)^/^ 

the total angular momentum of the binary reads 



(7) 
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Table 1. New physical parameters of some contact binaries. 
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6 


1.212 


0.852 


1.39 


0.70 


0.405 


7 


1.70 


0.83 


3.33 


0.71 


0.166 


8 


1.29 


0.95 


1.75 


1.05 


0.482 


9 


1.20 


0.71 


1.97 


0.66 


0.29 


10 


1.67 


0.61 


2.96 


0.50 


0.082 


10 


1.369 


0.797 


1.83 


0.71 


0.302 


11 


1.36 


0.62 


2.62 


0.54 


0.610 


12 


1.570 


1.220 


2.320 


1.44 


0.566 


13 


1.23 


0.75 


1.19 


0.54 


0.324 


14 



References in Table 1: (1) Gazeas et al. 2005; (2) Zola et al. 2005; (3) 
Zhang & Zhang 2004; (7) Yakut et al. 2005; (8) Qian et al. 2005; (9) 

Zhang & Deng 2005; (13) Ozdemir et 



Gazeas et al. 2006; (4) Terrell et al. 2006;(5) Yang et al. 2005; (6) 
Lu et al. 2007, (10) Szalai et al. 2007; (11) Lu 1991; (12) Zhang, 
al. 2002; (14) Djurasvic et al. 2001 



J = 54.23M«/^ A^/Vl (1 - m) + (1 + ^i^f) I , (8) 
in which 



where is a mass fraction of the primary of the system, 
and ri,2 the relative ra dii of both components. According to 
lEggletonl (|l98l l2006l l. they can be approximately written 
as 



fi = Ml /M, 



(9) 
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ri 



0.49{(1-m)/m} 



-2/3 



0.6{(1 - + ln{l + {(1 - ^^)/|.y'^'} 



(10) 



r2 ^ / 1 — A* \ 0.46 



(11) 



According to iBinnendiikI l|l97(t ). the contact condition can 
be expressed as 



A = ?±±^ = 1.327?i{l + 



0.76 



(12) 



Mass-radius relation for the primary of W- and A-type 
systems can be expressed as 

(13) 



Ri = c-Mf = crtj! M 

Inserting the equations (10)-(13) into equation (8), one can 
obtain 



J = 62.31cJ''^f(^)M 



in which 



/(m) 



(2+/3)/2 



then 



lOgJ: 



3 + /3 



^{l + (i^)i.92 



logAf + c(/i). 



(14) 



(15) 



(16) 



Based on Eqs. (14) and (16) in lAwadaUa fc Hannal l|2005l ). 
/3 can be derived to be 0.62 for A- and W-type systems, 
and Cr can be derived to be 1.05 (for W-types) or 1.23 (for 
A-types). Then, Cr is taken to be 1.23 since A-types are 
located near the two boundaries oi M ~ J diagram, and p is 
taken to be 0.5 (g = 1.0) and 0.93 {q — 0.075), respectively. 
Equation (16) is also shown in Figure 3 as two dashed hues 
with a slope of 1.81. It is seen in Figure 3 that all observed 
W UMa systems are located in a strip limited by the two 
dashed lines and the observed systems with extreme mass 
ratios, such as V870 Ara, CK Boo, FP Boo, V776 Cas, SX 
Crv, FG Hya, DZ Psc, AW UMa, GR Vir are on or near the 
lowest boundary (corresponding to a mass rat io close to a 
cutof f mass ratio of 0.076 for W UMa systems in iLi fc Zhani3 
l2006l ). The mass ratio of W UMa systems is decreased from 
the upper boundary to the low er boundary. 

lAwadalla fc Hannal (|2005 !') have computed only the or- 
bital angular momentum (OAM) of contact binaries and 
plotted it on a M — J diagram, then fit a quadratic function 
to it. However, a linear function can fit the distribution of 
the total angular momentum vs mass in their logarithmic 
scale for the observed W UMa systems well in present work, 
and a least-squares solution results a relationship between 
the total angular momentum and mass as the following 



logJ = 2.06(8)logi\f + 1.10(2), 



(17) 



with correlation coefficient r = 0.91 and standard deviation 
sd = 0.10. Equation (17) is also plotted in Figure 3 as a 
solid line with a higher slope than those of the two dashed 
lines based on Equation (16). This suggests that the mass 
ratio might be decreased with the decrease of the total mass 
for the observed W UMa systems. 

The distribution of mass ratio vs the total mass for 
the above systems is shown in Figure 4. It is seen in Fig- 
ure 4 that the W UMa systems seem to be located in a 



OS 



• A— type 
o W— type 

* A— type with low q 
A W— type with low q 



j8=0.62 



0.2 0.4 
logJl/(M„) 
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Figure 3. Total angular momentum of W UMa systems based on 
equation (6) as a function of total mass. Dashed lines represent 
the theoretical angular momentum distributions [based on equa- 
tion (16)] of W UMa systems with mass ratio of 1.0 and 0.075, 
respectively. Solid line represents a fitting result of observations, 
i.e. equation (17). 




¥(MJ 



0'' 

Figure 4. Mass ratio of W UMa systems as a function of to- 
tal mass. Full and open dots represent the A-type and W-type 
systems, respectively. 



strip and there is a tendency for a decreasing total mass M 
with dec reasing ma s s rati o q, which is similar to the result 
given bv lvan't Veeij (|l996l ). Meanwhile, as seen from Figure 
4, the mass ratios of W-type systems are located in a re- 
gion from 0.3 to 0.7 and the mass ratios of A-type systems 
seem to be located in two separated regions {q < 0.5 and q > 
0.7). Meanwhile, it is seen in Figure 3 that there is no obvi- 
ous difference between the angular momentum distributions 
of W-types and A-types. This is in good agreement with a 
result that the observed spectral type does not indicate a 
differe nce in evolutionary state between A- and W-type sys- 
tems l|Van Hammelll983 ). and A- and W-type systems are 
probably in the different stages of TRO. 
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The difference in the slope of equations (16) and (17) 
can be attributed to many causes, such as rough assump- 
tions, mean values of coefficients taken from empirical mass- 
radius relations, observational selection effects, initial mass 
and period distributions, and all effects including the effect 
of dynamical evolution. The rough assumptions and the co- 
efficient Cj, might only change the location (up or down) of 
the theoretical distribution given by equation (16) in Figure 
3, rather than its slope. The slope of equation (16) is only 
influenced by the coefficient (5 of empirical mass-radius rela- 
tion. However, equation (16) with a higher slope than equa- 
tion (17) r equires /3 > 1.0 which is mu c h larger t h an a v alue 
derived bv lAwadalla fc Hannal (120051 ). iRucinskil l|2006l ) has 
found that the contact binaries discovered by ASAS lacks 
those with photometric amplitudes < 0.4 mag. This might 
imply that a telescope has a detection threshold and some 
W UMa systems with smaller mass ratios will be difficult to 
be detected because of the very sraall amplitude of the corre- 
sponding light curves (|van 't Veeij|l978l : lRucinskill200lh . Fur- 
thermore, spectroscopic mass ratios were always determined 
for the best observable double lined spectroscopic W UMa 
systems. This sel ection favours the high mass ratio binaries 
l|van't Veedll978l '). Figures 3 and 4 show that the W UMa 
systems with smaller total masses usually have smaller mass 
ratios. Therefore, the observational selection effects possibly 
lead the slope of equation (17) to be decreased, ft is believed 
that the W UMa systems are formed fro m deta, c hed b ina- 
ries through AML due to MSW or through iKozail (Il963 ') cy- 



W UMa systems. In fact, it is consist ent with not only the 
prediction of the theoretical mode l s (ILucv 19761: Flanneryl 
19761 : iRobertson fc EggletonI Il977l : iLi. Han fc Zhand l2004l . 



20051 ) b ut also the observ ation al mass ratio d istributions 



cle (|Eggleton fc Kiseleval200ll : |Paczvriski et al.l2006l '). Then, 
the initial mass and period distribution would affect the 
time dependent formation process of these systems (i.e. the 
lifetime in the pre-contact stage). The lifetime in the pre- 
contact would affect mass, mass ratio, and nuclear evolution- 
ary degree of W UMa contact binaries just formed from the 
detached ones if their progenitors suffer mass loss and AML 
owing to MSW. The W UMa systems which originate from 
detached or semi-detached ones with main sequence compo- 
nents and short periods, usually spend a short t ime in pre- 
contact and have a high mass ratio of about 0.7 l|Mochnackil 
ll98ll : lLi. Han fc Zh ang 2004, 2005). W UMa contact bina- 
ries with a relatively long lifetime in pre-contact are usually 
formed from the detached ones which have a long period 
and a small radius ratio of the components to their Roche 
lobes. Since the more massive component evolves faster than 
the less massive one, the radius ratio of the secondary to 
its Roche lobe remains a relatively low value when the pri- 
mary begins to fill its Roche lobe so that the more mass 
is required to transfer from the primary to the secondary 
in the duration from the semidetached binary to a W UMa 
contact binary, with result that W UMa contact binaries 
should have a higher mass ratios. In addition, a lower rela- 
tive energy transfer rate is req uired for W UMa sy stems with 
a higher evolutionary degree ()Liu fc Yandl2OO0l i. the lower 
energy transfer rate should cause a smaller expansion of the 
radius of the secondary. In this case only the systems with 
a higher mass ratio might satisfy the requirement of Roche 
geometry. Therefore, W UMa systems with a relatively long 
lifetime in pre-contact should have a higher mass ratio and 
less mass. This would also lead the slope of equation (17) to 
be decreased. 

A higher slope of equation (17) compared with equation 
(16) is probably caused by the dynamical evolution of the 



given by Ivan 't Veeil i 19781 ) and iRucinskH (|200ll ) , who have 
argued that the W UMa systems with low mass ratios are 
much more than those with high mass ratios. If this is true, 
the W UMa systems would merge into the fast-rotating stars 
within 10^ — 10* yrs when th eir mass ratios is d ecrea sed to 
the c utoff mass ratio (0.076 iLi fc Zhan j |2006| . 0.09 iRasiol 
Il995h. Furthe r more, there is another possibility as suggested 
bv IPaczviiskil |2007|) that W UMa type binaries, especially 
for those with extreme mass ratios have secondaries in an 
advanced evolutionary stage as Algol-type ones. However, 
the secondary with a much larger radius than it would be 
on the main sequence is not only a result of its advanced 
nuclear evolution but also a result of energy transfer from 
the primary to it. If so, W UMa systems with extreme mass 
ratios and those with high mass ratios should have the dif- 
ferent formation processes. 



4 DISCUSSIONS AND CONCLUSIONS 

The status of the primaries of W UMa systems (including 
A- and W-type systems) predicted hj M — R diagram seems 
to be more advanced than that indicated by M — L diagram. 
Meanwhile, the primaries are almost located on ZAMS line 
for A- and W-type systems in Figures lb and 2b. The ap- 
pearance might be caused by the energy transfer from the 
primary to the secondary in W UMa systems, and the lu- 
minosity of the primaries of W UMa systems is affected by 
energy loss more greatly than their radii. In addition, the 
secondaries are significantly deviated from ZAMS stars in 
M — R and M — L diagrams. This might also be the result 
of a large amount of energy transfer from the more massive 
to the less massive component in W UMa systems. 

Since a large amount of energy is transferred from the 
primary to the secondary in W UMa systems, both com- 
ponents of W UMa systems can not achieve thermal equi- 
librium. At the beginning of contact stage, the efficiency of 
energy transfer is extremely low since the common envelope 
is too thin to undertake so much energy transfer, with the 
result that the mass is still transferred from the primary 
to the secondary as in semidetached stage. The decrease of 
orbital separation due to mass transfer should lead to the 
continued increase of the contact depth and energy trans- 
fer. The added mass on to the secondary would stop when 
the significant energy is transferred from the primary to the 
secondary. This followed by mass transfer back to the pri- 
mary until the decrease in radius of the secondary caused 
by mass loss can not compensated by its increase caused by 
energy transfer. Once this has occurred, the full efficiency 
of energy transfer is lost. Then, the secondary contracts and 
the contact depth decreases rapidly although the mass is 
transferred from the primary to the secondary in this stage. 
The secondary breaks contact and continues to collapse to- 
wards a main sequence equilibrium configuration, with its 
temperature and luminosity falling rapidly. The luminosity 
of the primary decreases during the semidetached phase, 
since the Roche lobe again contracts to prevent its free ex- 
pansion towards thermal equilibrium. Since the process of 
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raising matter up through the star, for transfer to the sec- 
ondary requires significant quantities of energy, at the ex- 
pense of the surface luminosity. Meanwhile, the timescales of 
the primary and the secondary are very different, because of 
mass accretion of the secondary its radius will increase be- 
fore it collapses to thermal equilibrium configuration. There- 
fore, the two components of W UMa systems are unlikely to 
be in thermal equilibrium, but each star attempts to reach 
thermal equilibrium. The attempt of the W UMa systems 
to reach a non-existent thermal equilibrium, coupled with 
Roche geometry, would force W UMa systems to undergo 
TRO on a period of abo ut thermal timescale of the primary 
(|Li. Han fc ZhaiiS |2004| . and references therein). Although 
the primaries of W UMa systems do not reach thermal equi- 
librium, they are near to the thermal equilibrium and exhibit 
the property that the luminosity is affected by energy trans- 
fer more greatly than its radius. But, the thermal timescale 
of the secondaries is much longer than the cyclic period so 
that the secondaries are evidently deviated from thermal 
equilibrium. 

Although Figures la and 2a show that the stage of evo- 
lution of A-type systems see ms to be higher than that of 
W-type ones as suggested by ILucv fc WilsonI (|l979l ') . How- 
ever, these objects are still within the main sequence band 
and we can hardly can speak of "evolved" obj ects, and this 
kind of binaries should last a very long time (|Van Hammd 
1 19821 ). Based on a simple assumption that the primaries 
of A- and W-type systems a re ZAMS stars, according to 
ICsizmadia fc KlagvivikI (|2004l ). the nuclear luminosity ratio 
of the secondary to the primary can be approximately ex- 
pressed as 



L2 ,M2s4,.e 4.6 



(18) 



Since the two components of W UMa systems have almost 
equal effective temperatures, i.e. Ti —T2, the observed lu- 
minosity ratio of them reads 

-^/2, observed 



L2 + AL _ Ra.2 



(19) 



where AL is the energy transfer rate between the compo- 
nents. Insert equation (19) into equation (20), one can obtain 
a relative energy transfer rate as 



AL 

"l7 



the relation of the relative energy transfer rate vs mass ratio 
is shown in Figure 5. It is seen in Figure 5 that the relative 
energy transfer rate of W UMa systems peaks at a mass ra- 
tio of about 0.55 which is just located in the region of mass 
ratio distribution of W-types, whilst the mass ratios of most 
A-types is significantly deviate from the maximum energy 
transfer rate of W UMa systems. Therefore, if the energy 
transfer is taken in account, we hardly can draw such a con- 
clusion that A-types are evolved more greatly than W-types 
so that A- and W-type systems have the the different evolu- 
tio nary status. This co nfirms the suggestions originally made 
bv lVan Hamiii^ l| 19821 ). which there is no a difference in evo- 
lutionary state between A- and W-type systems. Combining 
with the different mass ratio distributions of A- and W-type 
systems, we conclude that A- and W-type systems might be 
in the different stages of TRO, i.e. A-type systems might be 
in the stages form contact evolution to semi-detached one or 




Figure 5. Relative energy transfer rate between both components 
as a function of mass ratio of W UMa systems. 



from semi-detached configuration to contact one. This is ex- 
actly consistent with a relatively high temperature difference 
between the two components in A-type systems. However, 
no a TRO model can predict the observed number ratio of 
A-type systems with lower mass ratios to those with higher 
mass ratios. This suggests that present treatments of the 
energy transfer might be inadequate. 

The distribution of angular momentum and mass ra- 
tio suggests that the smaller is the total mass of the ob- 
served W UMa systems, the smaller is their mass ratios. It 
might be caused by the dynamical evolution of the W UMa 
systems, or by mass ratio re versal of the progenitors of W 
UMa systems (IStepieh| 2006h. espec ially for those with ex- 
treme mass ratios I Paczviiskil I2OO7I ). In addition, the first 
pos sibili ty is c onsi stent with the present theoretical mod- 
els (lLucv]il976i:iFlannery 1976 : Robertson fc Eg glcton 197^; 
iLi. Han fc Zhand 12004020051) and the observed mass ratio 
distribution (jvan't Veejll974 iRucinskil [ioOlh . i.e. the dy- 
namical evolution would lead W UMa systems to evolve into 
those with smaller mass ratios, then lead them to merge into 
fast-rotating stars owing to tidal instability. However, the 
second possibility might exist and it indicates that W UMa 
systems with extreme mass ratios and those with high mass 
ratios might have different formation processes. Meanwhile, 
the relative importance of these two possibilities requires 
further study through more observations. 

At present, it is not clear whether or not there is an evo- 
lutionary link between A- and W- type systems. If there is 
an evolutionary link between A- and W-type systems, based 
on the evolutionary tendency of mass ratio of W UMa sys- 
tems, we can conclude that A-type systems with relatively 
high mass ratios would probably evolve into W-type sys- 
tems, then evolve into A-type ones with relatively low mass 
ratios, and the opposite direction seems to be impossible, 
because it requires some W UMa systems with mass ra- 
tios close to 1 and different temperature between the two 
components. However, up to now there are in total of four 
contact binaries which are found to have mass r atios close 
to 1. they are V701 Sco feell fc Malcoiinl Il987l). CT Ta u 



(|Plewa fc Wlodarczvklll993l). V803 A ql (|Samec et al.lll993l ). 
and WZ And (jZhang fc Zhanall2006l ). Moreover, only a sys- 
tem WZ And in which the components with almost equal 
mass have different temperature between the components is 
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found, it might be unequal mass system through mass trans- 
fer and mass reverse. In each of the other three systems, the 
components with equal mass have almost identical temper- 
ature between the components. They could be regarded as 
congenital twins, rather than an outcome due to mass trans- 
fer during the evolution. 
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